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ABSTRACT 


A  technique  known  as  the  Taylor  Series  Correction  Method  (TSCM)  for 
extracting  the  stress  intensity  factor  from  photoelastic  data  is  reviewed. 
The  need  for  "artificial"  flaws  is  identified  and  an  approach  due  to 
Savin  is  used  to  evaluate  the  near  field  effects  of  various  practical 
flaw  shapes  upon  the  apparent  stress  intensity  factor.  Using  the  Sneddon- 
Srivastav solution  for  a  line  crack  in  a  finite  width  plate,  the  constric¬ 
tion  of  the  singular  zone  is  demonstrated  as  the  crack  tip  approaches  the 
free  edge.  Results  indicate  that  care  must  be  taken  in  applying  TSCM  to 
obtain  photoelastic  data  at  appropriate  distances  from  the  crack  tip. 


LIST  OF  SYMBOLS 


Polar  coordinates  as  defined  in  Figure  1. 

Cartesian  coordinates  as  defined  in  Figure  1. 
Nondimensional  coordinates  x/w,  y/w  as  shown  in 
Figure  9. 

Crack  length  or  depth  (for  surface  flaws). 

Width  of  symmetrical  notch  at  center  (Figure  4). 
Specimen  thickness. 

Width  of  finite  width  strip. 

Notch  angle  (Figure  4)  or  elliptical  flaw  angle. 
Complex  variable  x  +  iy  (Figure  A-l). 

Complex  variable  7  +  iy  (Figure  9). 

Radius  of  curvature  of  notch  at  tip. 

Fringe  order. 

Number  of  data  points. 

Material  fringe  value. 

Maximum  shear  stress  in  plane  perpendicular  to  the 
crack  border. 

Remote  tensile  stress. 

Regular  or  non-singular  stress  parallel  to  the  crack 
at  the  crack  tip. 


Kj,  SIF 
KAP 


Jq» 

V  Ti 

G(t) 


«o(c) 

$1 »  4*2 

4>»  'P 

A. 

i 


TSCM 

a 


Mode  I  stress  intensity  factor. 

Apparent  stress  intensity  factor  along 

IIiqX 

0  =  |  in  Figure  1 . 

Stress  components  in  a  plane  perpendicular  to  the 
crack  border. 

Stress  components  nondimensional ized  with  respect  to 
the  remote  tension  a. 

Bessel  functions  of  the  first  kind. 

Modified  Bessel  functions  of  the  first  kind. 

Solution  of  Fredholm  integral  equation  (6). 

Complex  variable  £  +  iq  defined  in  Figure  A-l . 

Mapping  function  which  maps  £  into  z  in  Figure  A-l. 
Mapping  angles  defined  in  Figure  A-l. 

Muskhelishvili  potentials. 

Least  squares  coefficients  for  G(t)  in  equation  (11). 
Coefficients  in  Taylor  Series  Correction  Method  of 
equation  (5). 

Taylor  Series  Correction  Method. 

Dummy  variable. 
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INTRODUCTION 


Although  the  use  of  photoelasticity  for  determining  stress  intensity 
factors  was  suggested  by  Irwin  [1]  as  early  as  1952,  extensive  employment  of 
the  method  has  been  slow  to  develop.  There  are  a  number  of  reasons  for  this. 
Results  of  Fessler  and  Mansell  [2]  and  Marl  off  et  al  [3]  revealed  the  need  for 
highly  accurate  measurements,  and  their  studies  were  reinforced  by  those  of 
Kobayashi  [4]-[6]  and  his  associates  in  their  work  with  propagating  cracks. 
Liebowitz,  Vanderveldt  and  Sanford  [7]  also  noted  certain  near  field  effects 
which  affected  stress  intensity  factor  (SIF)  determination.  In  order  to  obtain 
valid  data  for  SIF  determination  one  must  take  the  data  in  a  zone  dominated 
by  the  singular  stresses.  This  zone  is  bounded  on  the  outside  by  a  zone 
strongly  affected  by  the  remote  boundaries  and  loads,  and  on  the  inside  by 
a  zone  dominated  by  very  local  effects.  These  very  local  effects  for  a  stress 
freezing  photoelastic  material  consist  of  two  types: 

1)  Non-linear  effects  which  accrue  from  finite  rotations  and  deforma¬ 
tions  near  initially  sharp  crack  tips. 

2)  Notch  geometry  effects  when  the  crack  is  simulated  with  a  notch, 

or  when  the  crack  opens  a  substantial  amount,  producing  blunting.  Beginning 
in  1969,  the  senior  author  and  his  associates  carried  out  a  series  of  stress 
freezing  photoelastic  investigations  [8]-[19].  These  studies  revealed  that, 
in  a  nunber  of  three  dimensional  problems  [13]-[15]  the  zone  dominated  by  the 
singular  stresses  was  severely  constricted  and  that  it  was  necessary  to  use 
data  outside  the  singular  zone  for  SIF  determination.  In  order  to  account 
for  remote  effects  included  in  this  way,  a  technique  employing  a  Taylor 
series  expansion  of  the  maximum  in-plane  shear  stress  to  account  for  the 
regular  part  of  the  stress  field  was  employed.  The  method  is  referred  to 
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as  the  Taylor  Series  Correction  Method  (TSCM).  Details  of  the  method  have 
been  described  in  [14]  and  [17]  and  its  applicability  to  three  dimensional 
problems  was  recently  discussed  in  [20].  (This  method  will  be  featured  in 
a  chapter  in  the  second  volume  of  the  SESA  Monograph  on  Experimental  Frac¬ 
ture  Mechanics).  Although  the  use  of  TSCM  has  extended  the  applicability  of 
photoelastic  SIF  determination  to  three  dimensional  problems  not  previously 
solved,  practical  difficulties  remain  in  fabricating  appropriate  crack  geome¬ 
tries,  producing  effects  upon  the  inner  zone  previously  mentioned.  More¬ 
over,  analytical  verification  of  the  constriction  of  the  singular  zone  due 
to  free  boundaries  or  three  dimensional  effects  has  not  been  demonstrated. 

In  a  very  recent  analytical  study  [21]  McGowan  and  Smith  have  shown  that  the 
non-linear  effects  in  the  inner  zone  are  much  too  near  the  crack  or  notch 
tip  to  be  picked  up  in  the  photoelastic  data.  However,  in  problems  where 
a  specific  crack  shape  is  desired,  an  artificial  crack  may  be  sawed  in, 
producing  a  relatively  large  root  radius,  or  the  photoelastic  material  in 
the  rubbery  range  may  stretch  enough  locally  near  the  crack  tip  under  even 
small  loads  to  produce  stress  field  disturbances  which  may  extend  into  the 
photoelastic  data  zone.  The  present  study  is  concerned  with  a  quantitative 
analytical  assessment  of  effects  of  this  type,  and  also  with  a  two  dimensional 
analytical  study  of  the  influence  of  free  boundaries  upon  the  singular 
zone- 


BASIC  THEORY 


The  Irwin  Equations  [1]  for  the  elastic  stress  field  near  the  tip  of 
a  plane  crack  are: 


*1 

(2rrT 


■1/2 


cos  {1  -  sin  sin  -  cl. 


=  (5r)1/2  005  *  {1  +  Sl"  ?  Sl" 


(1) 


rxy 


sin  |  cos  |  cos 


with  the  notation  given  in  Figure  1.  This  discussion  will  be  confined  to 
plane  problems  so  that  complete  solutions  can  be  employed  in  all  mathema¬ 
tical  models.  Application  of  the  method  to  three  dimensional  problems  is 
described  elsewhere,  [14],  [19],  [20].  The  quantity  measured  photoelasti- 
cally  i‘s  the  fringe  order  which  is  proportional  to  the  maximum  in-plane 
shearing  stress.  The  local  value  of  this  latter  quantity  may  be  computed 
by  substituting  Equations  (1)  into: 

Tmax  =  ?  -oy  )c  +  4^  (2) 

Then  evaluating  rmax  along  the  convenient  direction  e  =  #/2  and  combining 
this  result  with  the  stress  optic  law 

T  =  Ilf. 

max  2t  (3) 

one  can  obtain: 


Kj  =  Ur)1/2 


(4) 


Equations  (1)  constitute  a  two  parameter  (Kj,ct0)  set  of  field  equations  and 
are  valid,  according  to  Irwin  in  the  singular  zone. 

When  one  uses  data  outside  this  zone,  the  following  equation  may  be 


used: 


max 


K  M  N 

_I_  l  /?  +  £  BMrN 
(8irr) ' ' 2  N=0  N 


(5) 


instead  of  Equation  (4)  where  the  first  term  arises  from  the  singular 
stresses  and  the  Taylor  Series  accounts  for  the  regular  part  of  the  stress 
field  including  aQ.  Use  of  Equation  (5)  is  made  in  the  TSCM  and  the  number 


of  degrees  of  freedom  of  the  system  can  be  greater  than  two  (i.e.Kj,  BQ,  B-j , 

&2*  •  •  •  Bn)  but  normal does  n°t  exceed  five  or  six. 

Although  one  normally  would  expect  to  plot  stress  (i.e.x  )  versus 

max 

distance  from  the  crack  tip,  the  authors  have  found  this  type  of  plot  to  be 
relatively  insensitive  to  changes  in  the  stress  intensity  factor.  Instead, 
an  apparent  stress  intensity  factor  is  defined  as: 

KAP  ■  <8->V2  («) 

and  a  normalized  value  of  this  quantity  is  plotted  against  (r/a)^. 
Extrapolation  of  this  curve  to  r/a  =  0  yields  the  experimental  SIF. 

NOTCH  EFFECTS 


Figure  2  is  a  microphoto  of  a  natural  flaw  under  a  field  of  remote 
tension  perpendicular  to  the  crack  plane  showing  how  a  "hairline"  crack 
opens  under  tension  into  a  geometry  approximating  a  narrow  ellipse.  This 
blunts  the  crack  tip  into  a  root  radius  of  about  5  x  10"3  mm.  An  assess¬ 
ment  of  the  influence  of  this  blunting  was  made  by  the  authors  in  an 
earlier  paper  [12]  by  computing  the  apparent  SIF  from  the  Ingl is-Kolosoff 
linearly  elastic  solution  for  an  elliptical  hole  in  an  infinite  plate 
under  tension  and  normalizing  it  with  respect  to  the  SIF  for  a  line  crack 
in  an  infinite  plate.  The  results  are  shown  in  Figure  3  and  reveal  a 
blunting  effect  which  is  negligible  only  if  data  are  taken  at  least  five  to  ten 
root  radii  from  the  crack  tip. 

Artificial  notches  may  be  made  in  a  variety  of  geometries.  A  series 
of  those  found  convenient  by  the  authors  are  pictured  in  Figure  4  together 
with  the  elliptical  notch  just  discussed.  Approximate  mapping  functions 
have  been  provided  by  Savin  [22]  which  lead  to  the  determination  of  the 
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Muskhelishvili  complex  potentials  [23]  from  which  complete  solutions  to 
these  geometries  may  be  obtained.  Such  procedures  are  detailed  in  Appendix 
A  for  the  geometries  (other  than  elliptical)  of  Figure  4. 

The  solution  for  a  rectangular  notch  was  obtained  as  a  special  case  of 
the  general  mapping  solution  presented  in  Appendix  A.  The  solution  is  exact 
for  the  region  of  the  z  plane  corresponding  to  the  unit  circle  in  the  x, 
plane.  This  region  is  not,  however,  an  exact  rectangle  but  has  curvilinear 
sides  and  corners  of  some  small  radius  of  curvature.  By  increasing  the 
number  of  terms  in  the  mapping  function,  rectangles  having  straight  edges 
and  corners  of  any  degree  of  sharpness  can  be  obtained.  The  limiting  case 
has  stress  singularities  at  the  corners. 

Figure  5  is  a  plot  of  the  normalized  apparent  SIF  for  rectangular 
notches  corresponding  to  Figure  3  for  elliptical  notches.  The  corresponding 
curve  for  the  limiting  case  of  a  line  crack  is  again  a  straight  line  but 
in  this  case  has  a  positive  slope  of  unity  which  results  from  the  fact  that 
CTq  =  a  for  uniaxial  tension.  It  should  be  noted  that  locally  the  two  solutions 
have  an  entirely  different  character  with  the  rectangular  notch  having  peak 
stresses  at  the  corners  and  the  ellipse  having  peak  stresses  at  the  notch  tip. 

A  comparison  of  Figures  3  and  5  shows  a  remarkable  similarity  between  the 
elliptical  notches  with  p/a  =  .01  and  .05  and  the  rectangular  notches  with 
b/a  =  .01  and  .05.  Specifically,  it  can  be  seen  that  the  notch  curves  begin 
to  coincide  with  the  line  crack  solutions  at  the  same  value  of  r/a  whether  the 
notch  is  elliptical  or  rectangular  if  we  associate  the  value  of  b/a  for  the 
rectangular  notch  with  that  of  p/a  for  the  elliptical  notch.  This  implies  that 
it  is  not  the  minimum  radius  of  curvature  which  is  important  but  the  relative 
width  of  the  notch  at  the  tip  (note  that  the  rectangular  notch  has  extremely 
small  radii  of  curvature  at  the  corners  and  in  the  limit  predicts  stress 
singularities) . 
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As  pointed  out  earlier,  the  number  of  terms  retained  in  the  mapping 
function  affects  the  radius  of  curvature  at  the  corners.  Figure  5  represents 
10  terms  in  the  mapping  function.  Figure  6  shows  the  effect  of  varying  the 
number  of  terms  on  one  of  the  curves  of  Figure  5.  It  is  obvious  that  unless  a 
very  small  number  of  terms  is  being  used  this  effect  is  secondary  on  the 
comparison  of  interest. 

The  case  of  angular  notches  was  also  studied  as  a  special  case  of  the 
mapping  solution  of  Appendix  A.  The  character  of  the  solution  is  somewhat  different 
as  illustrated  by  the  sharp  tip  (p  <  5  x  10"3mm)  rhombus  solution  shown  in 
Figure  7.  Even  for  very  shallow  angles  the  notch  solution  is  seen  to  be 
somewhat  below  the  line  crack  solution.  At  the  origin  the  curves  converge  which 
indicates  that  the  stress  intensity  factor  is  the  same  as  that  for  the  line  crack. 
This  is  further  substantiated  by  the  work  of  Gross  and  Mendelson  [24].  How¬ 
ever,  away  from  the  origin  where  experimental  data  would  be  gathered  the  solutions 
differ  by  about  10%  for  a  30°  notch.  Hence  a  linear  extrapolation  would  tend  to 
underestimate  the  stress  intensity  factor.  This  effect  is  partially  due  to 
the  large  value  of  b/a  involved  in  the  rhombus. 

The  same  comparison  for  a  blunted  triangular  notch  is  shown  in  Figure  8. 

This  type  of  notch  is  easily  obtained  in  practice  by  grinding  the  corners  from 
a  standard  circular  saw  blade. 

This  modification  is  very  useful  since  the  solution  behaves  essentially 
like  a  rectangular  notch  of  width  2p  rather  than  2b  (compare  Figures  5  and  8). 

As  in  the  previous  solution  (Figure  7)  the  notch  stresses  are  somewhat  below  the 
crack  stresses;  however,  for  b/a  <  .10  the  difference  is  quite  acceptable. 

FREE  BOUNDARY  EFFECTS 

The  effect  of  a  free  boundary  on  the  stress  field  local  to  a  crack 
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tip  has  been  of  considerable  interest  in  the  field  of  fracture  mechanics 
for  some  years.  Numerous  analytical  and  numerical  investigations  have  been 
carried  out  with  the  specific  purpose  of  determining  free  surface  magnifi¬ 
cation  factors  for  cracks  approaching  a  free  boundary.  However  since  the 
stress  intensity  factor  can  be  extracted  from  these  solutions  with  much 
less  labor  than  a  complete  solution  requires,  analytical  solutions  are 
limited  almost  invariably  to  Kj  itself  and  give  no  information  about  higher 
order  terms  in  the  stress  components.  This  causes  two  significant  problems 
in  attempting  to  verify  a  crack  solution  experimentally.  First,  in  the 
absence  of  the  complete  solution,  it  is  impossible  to  predict  a  priori  the 
size  of  the  zone  dominated  by  Kj  and  <Jq  for  use  in  the  two  parameter  tech¬ 
niques  previously  discussed.  Second,  it  is  very  difficult  to  judge  whether 
experimentally  observed  trends  are  valid  without  a  knowledge  of  the  quali¬ 
tative  character  of  the  analytical  solution  away  from  the  crack  tip. 

Most  numerical  techniques  are  not  well  suited  for  studying  this  effect 
since  accurate  results  are  desired  in  the  neighborhood  of  both  a  singulari¬ 
ty  and  a  free  boundary.  Therefore  any  technique  which  does  not  make  special 
provision  for  the  singularity  or  which  collocates  on  the  boundary  is  unde¬ 
sirable. 

A  convenient  geometry  for  studying  this  effect  is  the  through  crack  in 
a  finite  width  strip  shown  in  Figure  9.  Although  this  problem  is  two  dimen¬ 
sional  the  qualitative  nature  of  the  solution  is  expected  to  model  very 
closely  the  behavior  of  a  three  dimensional  surface  flaw  approaching  a  free 
boundary  since  the  local  behavior  of  the  three  dimensional  flaw  at  maximum  depth 
is  close  to  plane  strain.  Sneddon  [25]  has  reduced  the  solution  of  this 
problem  to  the  solution  of  the  following  Fredholm  integral  equation  of  the 
second  kind. 
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(6) 


G(t)  =  t  +  t 


fa 

'o 


G(u)H(u,t)du 


where  a  is  the  half  crack  length  a  nondimensional ized  with  respect  to  the 
half  width  of  the  strip  w  and  the  kernel  is  given  by: 


H(u,t)  = 

where 


roo 


xe"x 

sinh(x) 


IQ(xu)  I 0 ( xt)  + 


2xh(u,x)h(t,x)  1 
2x  +  sinh(2x)  J 


dx 


(7) 


h(u,x)  =  [x  coth(x)  -  1]  Iq(ux)  -  uxl^(ux) 

and  Iq  and  I-j  are  modified  Bessel  functions  of  the  first  kind.  This  equa¬ 
tion  has  been  solved  numerically  by  the  method  of  successive  substitutions 
utilizing  a  subroutine  recommended  by  Squire  [26]  for  the  kernel  evaluation. 
The  subroutine  is  designed  for  infinite  integrals  and  gives  two  independent 
estimates  of  the  integral  based  on  a  20  point  Gaussian  formula  and  a  41  point 
Kronrod  rule.  The  difference  between  the  estimates  was  of  the  order  of  10”^ 
for  this  particular  integrand  and  can  be  taken  to  be  a  conservative  estimate 
of  the  error  in  the  integrand.  The  technique  for  solving  (6)  consisted  of 
splitting  the  t  axis  into  N  discrete  points  t.,  letting  Gn(t.)  =  t.  and 
determining  successively  G.|,  Gg»  .  .  •  G^  until  convergence  was  achieved  by 
use  of  the  formula 


G 


i+1 


‘j  *  ‘j 


fa 

G.(u)  H(u,t.)  du 
Jq  1  J 


(8) 


The  definite  integral  was  evaluated  by  a  subroutine  from  [26]  very  similar 
to  the  one  previously  discussed  only  for  a  finite  interval.  The  solution 
was  found  to  be  relatively  insensitive  to  the  number  of  divisions  along  the 
t  axis.  The  results  for  G(t)  are  plotted  in  Figure  10  for  various  values 
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_  1  /?  _ 

of  a/w.  The  stress  intensity  factor  is  given  by  Kj/o  (ira)  '  =  G(a)/a  as 
derived  from  reference  [25].  It  was  computed  for  Values  of  a/w  ranging 
from  0.1  to  0.9  and  agreed  within  1/2%  with  Sneddon's  values. 

Having  determined  G(t)  the  stresses  are  obtained  by  superimposing  the 
following  two  stress  systems  as  reduced  from  reference  [25]. 

Stress  System  I 


a  +  a 

_X _ £  = 

2 


<j>(p)e“py  cos(px)  dp 


a  -  a 

■Vs 


Jo 


p4>Cp)emf>y  cos(px)  dp 


(9) 


t  =  y 
xy  J 


JO 


where  <f>(p)  =  -p 


p4>(p)e"py  sin(px)  dp 


G(u)  Jq(pu)  du,  Jq  is  the  zeroeth  order  Bessel 


! function  of  the  first  kind,  the  stresses  have  been  nondlmensionalized  with 

i 

respect  to  the  remote  stress  a,  and  x,  y  are  nondimensional  coordinates  as 
shown  in  Figure  9. 

Stress  System  II 


*00 

g(p)  cosh  (px)  cos(py)  dp 

Jo 


2 


r°° 

[h(p)  cosh(px)  +  px  g(p)  sinh(px)]  cos(py)  dp 


JO 

*00 

[h(p)  slnh(px)  +  px  g(p)  cosh(px)]  sln(py)  dp 

-0 


(10) 
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where 


g(p) 

h(P) 


jrr^STOT[2C2p  +  C1(3-20t.-^)] 

2p  +  sinh(2PT  [  2C1 (1  -  2p  +  p2  +  e  Zp>  +  C2P  (1  -  2p  +  e  2p)  ] 


C-|  (p) 
C2(p) 


G(u)  In(pu)du 
J0  u 


a 

G(u)ul-,  (pu)du 

0 


and  IQ  and  I-j  are  as  previously  defined. 

Stress  system  I  is  singular  and  for  the  special  case  G(u)  =  u  reduces  to 

the  solution  for  a  through  crack  in  an  infinite  body  in  tension.  For  this 

case  (|>(p)  =  -a  J^pa)  which  indicates  that  the  integrand  in  (9)  can  be  expected 
to  be  highly  oscillatory.  In  addition,  since  e“py  decays  very  slowly  for 
the  region  of  interest  (small  y),  it  is  apparent  that  a  direct  numerical 
approach  to  this  system  of  equations  will  be  unsatisfactory.  Instead  an 
expression  of  the  form 

6<t>  ■  I  A/"'1  (ID 

n=l 

was  fitted  to  the  numerical  values  of  G(t)  by  the  method  of  least  squares. 

The  use  of  equation  (11)  allowed  the  evaluation  of  equations  (9)  in  closed 

form  as  shown  in  Appendix  B.  It  also  permitted  the  closed  form  evaluation  of 

C-j  and  C2  after  which  stress  system  II  could  be  evaluated  using  standard 
techniques.  The  integrand  of  stress  system  II  was  also  oscillatory  due  to 
sine  and  cosine  terms,  but  it  converged  quite  rapidly  due  to  the  term  sinh(2p) 


in  the  denominator. 
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Since  the  solution  for  G(t)  (Figure  10)  is  nearly  linear  even  for  large 

values  of  a/w  it  was  found  that  three  terms  in  (11)  were  sufficient.  The 

approximation  of  G(t)  had  no  effect  on  the  free  edge  boundary  conditions 

7  (+  l,yl  =  x  (+  l,y)  =0  since  they  follow  from  the  form  of  g  and  h  in 
x  xy 

equation  (10).  The  shear  stress  along  both  planes  of  symmetry  also  vanished 
identically  so  that  the  only  effect  of  the  approximation  was  that  residual 
values  of  7  remained  on  the  crack  surface.  The  average  of  these  residuals 

y  -  n  -2  2 

varied  from  10”-  N/mnr  for  a/w  =  0.1  to  nearly  10  N/mm  for  a/w  =  0.9. 

In  all  cases  they  are  small  relative  to  the  applied  tension  which  is  unity. 

Figure  11  is  a  nondimensional  plot  of  K^p  vs.  r^  based  on  numerical 

results  of  the  finite  width  strip  solution.  The  two  parameter  singular  zone 

is  the  linear  portion  of  each  curve  near  the  origin.  As  illustrated  by  the 

dashed  line  passing  through  the  point  of  tangency  drawn  from  the  origin 

along  each  curve,  the  singular  zone  shrinks  rapidly  as  the  crack  approaches 

the  free  edge  of  the  strip.  Perhaps  even  more  significant  is  the  slope  of 

the  curves  at  the  origin  which  is  proportional  to  <?q.  As  the  crack  ti|p 

I. 

approaches  the  free  boundary  the  slope  increases  significantly.  This  jis 

demonstrated  in  Figure  12  which  is  a  non-dimensional  plot  of  aQ  and  Kjj  vs. 

a/w.  It  demonstrates  that  as  a/w  approaches  unity  not  only  is  Kj  magnified 

but  also  <7q  and  that  the  increase  in  aQ  even  exceeds  that  in  Kj.  This  implies  that  aQ 

is  primarily  due  to-: the  interaction  between  the  singular  stress  field  and  the 

free  boundary  and  has  little  to  do  with  the  remote  stress  field  for  problems 

with  near  boundaries.  This  is  an  important  observation  for  most  surface  flaw 

problems  of  current  practical  interest  which  invariably  involve  free  boundaries. 

CONCLUSIONS 

From  this  study  the  following  general  conclusions  can  be  made: 
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A  -  Notch  Effects 


1)  When  possible,  natural  flaws  should  be  used  to  simulate  the  line 
crack.  In  such  cases  data  gathered  outside  of  r/a  =  lOp/a  will  ordinarily  be 
free  of  notch  effects. 

2)  Rectangular  notches  can  be  used  to  simulate  line  cracks  provided 
that  data  are  gathered  outside  of  some  minimum  radius.  This  radius  can  be 
judged  from  the  elliptical  notch  solution  (Figure  6)  provided  that  the  half 
width  of  the  rectangular  notch  is  associated  with  the  radius  of  curvature 
of  the  elliptical  notch. 

3)  Blunt  triangular  notches  can  be  used  to  simulate  the  line  crack 
provided  that  b/a  ~  .10  and  data  are  gathered  outside  of  a  specified  value  of 
r/a  which  can  be  determined  from  the  elliptical  notch  solution  by  associa¬ 
ting  the  half  width  of  the  triangular  notch  at  the  tip  with  the  radius  of 
curvature  of  the  elliptical  notch. 

4)  Sharp  angular  notches  will  tend  to  significantly  underestimate  the 
stress  intensity  factor  unless  data  are  taken  extremely  close  in. 

B  -  Free  Boundary  Effects 

5)  The  singular  zone  is  shown  to  shrink  rapidly  as  a  crack  approaches 
the  free  edge  of  a  strip. 

6)  When  a  crack  is  near  a  free  boundary,  Oq  depends  primarily  upon 
the  interaction  between  the  singular  stress  field  and  the  free  boundary,  and 
is  not  strongly  influenced  by  the  remote  stress. 

All  of  the  foregoing  conclusions  are  specifically  directed  towards  clari¬ 
fying  the  determination  of  SIF  values  by  means  of  the  TSCM  utilizing  stress 
freezing  photoelasticity  which  involves  no  plastic  zone.  However,  there  is 
a  very  local  non-linear  region  [21]  which  is  expected  to  produce  effects  in 
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the  singular  zone  which  would  be  similar  to  those  associated  with  small 
scale  yielding  in  structural  materials. 
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APPENDIX  A 
MAPPING  SOLUTION 

It  is  well  known  that  the  stresses  around  a  hole  in  a  homo¬ 
geneous  isotropic  plate  can  be  obtained  by  the  complex  variable  tech¬ 
nique  introduced  by  Muskhel ishvili  [23].  Numerous  problems  involving 
square,  rectangular  and  triangular  holes  have  been  solved  by  Savin 
[22].  For  the  most  part  these  solutions  contain  only  a  few  terms  in 
the  series  expansion  of  the  mapping  function  and  hence  involve  curvi¬ 
linear  sides  and  rounded  off  corners.  With  the  advent  of  the  digital 
computer  it  has  become  feasible  to  automate  the  process  of  determining 
the  Muskhel ishvil i  potentials  <j>  and  and  hence  to  arrive  at  solutions 
containing  arbitrarily  sharp  corners  and  straight  edges.  Thus  it  is 
possible  to  determine  the  effect  of  the  geometrical  approximations 
involved  in  truncating  the  series  mapping  function. 

For  the  present  purpose  of  determining  the  effect  of  various 
types  of  notches  on  the  estimation  of  the  mode  I  stress  intensity 
factor  it  is  desirable  to  consider  the  axisymmetric  eight  sided  polygon 
of  Figure  A-l .  Once  the  general  solution  is  obtained  for  this  region 
the  following  special  cases  can  be  evaluated: 

a)  Rectangular  notch  (a  =  0,  <f>}  =  <f>2  ?  0) 

b)  Rhombus  (a  f  0,  <j>-|  =  0,  <j>2  =  |-) 
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c)  Rectangular  notch  with  sharp  triangular  tip 
(a  f  0,  <f>]  =  0,  <(>2  t  0) 

d)  Rectangular  notch  with  blunted  triangular  tip 
(a  f  0,  <(>1  t  0,  <f>2  f  0) 

To  obtain  the  mapping  of  Figure  A-l ,  the  Schwartz  Christoffel  transfor¬ 
mation  is  first  applied  which  maps  the  exterior  of  the  polygon  in  the 
z  plane  into  the  exterior  of  the  unit  circle  in  the  ?  plane.  This 
transformation  is  given  by  (see  [22]) 


‘c 

[l-f] 

or|-l 

(i-  *r] 

a2~^ 

*  •  • 

r  n'l 

i  -  r 

f  ' 

w  * 

l 

dt 


+  complex  constant.  (A-l) 

Referring  to. Figure A-l  it  is  apparent  from  symmetry  considerations  that 
A'  =  -E'  =  e_i(f)l 

_i<f>2 

B'  =  -F1  =  e 

-  Up 
C'  =  -G'  =  F'  =  -e  c 

1*1 

D’  =  _h'  =  E1  =  -e  1 


and 


3-a 

a]  =  a4  =  a5  =  a8  =  2 

1+a 

a2  =  a3  -  a6  =  a?  -  2 


Hence  equation  (A-l)  takes  the  form 
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z  =  R 


1  - 


(1 -o)/2 


21  ♦! 


o/2  , 


1  - 


dt  +  complex  constant. 


e“2i<J>2i  (l-a)/2 

1H 


(A-2) 


Since  |e19|  =  1  and  |t|  >"1,  this  equation  can  be  approximated 
by  expanding  each  of  the  factors  in  a  binomial  series  and  multiplying 
the  resulting  polynomials.  Then  by  integrating  the  resulting  poly¬ 
nomial  and  by  inversion  (i.e.  by  replacing  t,  by  ^-)  the  desired  mapping 
is  obtained. 

As  a  result  of  symmetry  it  is  found  that  the  required  mapping 
is  of  the  form 


z  =  R 


U 

c 


N 

I 

n=l 


2n-l 


=  «(0 


(A-3) 


where  the  coefficients  R,  Cn  are  real  and  the  expansion  has  been 
truncated  after  the  Nth  term.  A  computer  program  was  written  to  carry 
out  the  above  calculations  for  given  values  of  4>i ,  4>2»  a  and  N. 

It  is  well  known  that  the  solution  of  the  first  boundary  value 
problem  of  elasticity  can  be  reduced  to  the  determination  of  two 
analytic  functions  <f>(z)  and  \p(z)  of  the  complex  variable  z.  The 
stresses  are  related  to  the  complex  potentials  <|>  and  <p  by  the  following 
formulae  where  primes  denote  differentiation: 


ax  +  ay  =  2[<f> '  (z)  +  <(>•  (z)] 

ay  -  ax  +  2ixxy  =  2[z<j>"(z)  +  V(z)] 


(A-4) 
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The  stress  boundary  conditions  using  this  formulation  take  the  form 

<j>(z)  +  z<J>'(z)  +  #(z)  =  f]  +  if2  on  C  (A-5) 

s 

where  f-j  +  i f2  =  i  (  [T-j(s)  +  ^(s)]  ds  is  the  force  transmitted 

s0J 

across  a  segment  of  the  boundary  (ds)  and  T-| ,  T2  are  components  of  the 
applied  tractions  on  the  boundary  C.  Under  the  transformation  z  =  u(?) 
(using  the  notation  <f>(6)  =  <J>(“(?))  etc.)  the  stress  transformation 
equations  yield  the  following  expressions  for  the  stress  components 
Op  *  Tp^,  cty  in  the  curvilinear  coordinate  system  corresponding  to 
the  contour  lines  p  =  constant  and  <j>  =  constant  where  z  =  wfpe1*^). 

Op  +  =  2[(J>(c)  +,<Kc)] 

2,2  _  (A-V 

<y<j>  -  ap  +  2ixp,  =  ~2§  —  Cw(c)  *'(0  +  (o'(c)«(c)] 

P  u  (?) 


The  stress  functions  <f>  and  can  be  separated  into  components 
<f>*>  1>*  representing  the  stress  state  in  the  absence  of  the  hole  and 
<J>o>  <Pq  due  to  the  presence  of  the  hole.  For  uniaxial  tension  in  the 
y  direction  <f>*  and  \p*  take  the  form 


.*  _  sL  _  o  /  , 

~  4  -  4  w(?) 
4>  =  2  =  2 


(A-7) 


The  potentials  4>o  an<i  can  be  expanded  as  follows  (where  only  odd 
powers  need  be  retained  due  to  symmetry  in  both  region  and  loading): 
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CO 


00 


*0  ■  £,  V= 

n=1 


+o  =  Bn? 
n=l 


(A-8) 


As  first  shown  by  Muskhelishvili  [23]  the  boundary  conditions  can  be 
written  in  terms  of  the  two  functional  equations: 


♦o(5)  +  2iri  u'(5)  *o'(«)  S-C  +  B0  2iri  Jp  0_ 5 
r 


f,  +  if2 


4*0  (?)  +  2iri 


__y(<jj)  .  ,  do_  _  J_  f  1  1  f2 

u'  (5)  ^0  a-£  2iri  0^  ^0 

r  Jr 


(A-9) 


where  r  represents  the  unit  circle.  For  uniaxial  tension  in  the  y 
direction  it  can  be  shown  that  (see  [22]): 


fo  “  0 


<1  *  Mi)  +  £(i)] 


(fl-10) 


Substituting  (A-3)  in  ( A— 1 0 )  and  noting  that  the  mapping  coefficients 


are  real  yields: 


rl  =  ~y 


1  .  ”  -2n-l  1  J  -  2n-l 

?  +  l  cn°  +  J  *  E,  cn0 
L  n-l  n=l 


(A-ll) 


But  o  =  a  ;  hence  (A-ll)  becomes 


f,  =  _  crR  1  +  ±  .  v  r  ^2n-l  -l-2nx 
rl  T  a  +  0  +  E  Cn(a  +  0  ) 

1_  n=l 


(A-12) 


Hence  the  right  hand  side  of  equation  (A-9)  becomes 
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1 

■  f]da  aR. 

N  2n-l 

2rri 

5-?  =  -  2 
r 

5  +  nf,  C"£ 

^  J 

(A-13) 


where  the  following  result  from  the  theory  of  residues  has  been  used: 


1 

2TTi 


(  onda 

— x - 

CT-C 


cn  n  >  0 
0  n  <  0 


(A-l 4 ) 


Utilizing  equation  (A— 3 )  and  noting  that  a )'(a)  f  0  and  hence 
>'(a)  t  0  for  single  valued  mapping  functions  we  obtain  the  identity: 


1 

a  + 


g(g)  = 

u'(ct)  _-j 


l  Cna 
n=l 


.2n-l 


N 


N-l 

=  Z  H„a 


2n-l 


+  Z  (2n-l )Cna 
n=l  n 


2-2n  n=l 


+  negative  powers  of  a 

where  the  H's  are  obtained  by  long  division  and  negative  powers 
of  a  are  unnecessary  because  they  do  not  contribute  to  the  Cauchy 
integrals.  Also  from  equation  (A-8)  we  have 


*0'(i)  =  £  (2n-l)  flnS2'2n  (A.,6) 

n=l  ' 

Multiplying  equations  (A-l 5)  and  (A-16)  yields 
w(a)  r  „  *2n-i 

=T(5)  <p0  ^  Kna  +  negative  powers  of  a  (A-17) 

where  and  A^  are  related  by  the  matrix  equation: 
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K-j  =  Bj j  Aj  (implied  summation)  (A-18) 

and  the  coefficient  matrix  is  given  by 

Bi j  =  (2j-l )  (no  summation).  (A-19) 

Making  use  again  of  relation  (A— 14 )  the  second  term  in  the  first  of 
equation  (A-9)  becomes 


N-!  9n  ,  ^  N-l  2n-l  N“1  -  J 

r  KnS2"-1  ^  -  I  K„c  *  f  W 

r  n=l  n=l  n=l 


(A-20) 


Substituting  equations  (A-8),  (A-13),  and  (A-20)  in  the  first  of 
(A-9)  yields  the  result  that  B0  =  0  and  provides  the  following  system 
of  equations  to  solve  for  the  An: 


A1  +  Blj  Aj  ’  '  2  +  Cl] 


A  +  B  •  A-  =  -  —  C 
Mn  Dnj  Ho  2  n 


n  =  2,  3,  ...  N-l 


(A- 21 ) 


Since  the  mapping  coefficients  Cn  are  real  for  the  cases  to  be  investi¬ 
gated  here  the  coefficients  An  will  obviously  be  real  also.  Hence 
equation  (A-21)  can  be  written  in  the  concise  form 


<5nj  +  V  Ao  ■  -  r  [sm  +  cn] 


(A-22) 


where  6^  is  the  Kronecker  delta. 
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Once  <f>0  is  known  ij>0  can  be  determined  from  the  second  of 
equations  (A-9)  which  can  be  reduced  to  the  following  form  (see 
Sokol nikoff  [28] »  equation  84.11): 


Ve)  '  2iri 


if. 


5  -  c 


6a 

0)  U 


N-l  _  l-2n 
*o’(t)  +  „=1  KnC 


(A-23) 


where  for  real  coefficients  Kp  =  Kn  and 


oj*(?)  =  R 


i+ 


N 

Z 

n=l 


*  o> (6 ) 


Substituting  equations  (A-3)»  (A-8),  and  (A-13)  in  (A-23)  we  obtain 
finally: 


N 


?  +  l 

n=l 


N-l 


+  z 
n=l 


KnC 


l-2n 


N 

?  +  z  Cn? 
n=l 


l-2n 


1 


N 

z  (2n-l )  Cn? 
n=l 


2n-2 


N  2n-2 

z  (2n-l )  Anc 
n=l 


(A-24) 


This  completes  the  formal  solution.  The  previous  calculations 
were  programmed  for  a  digital  computer  which  computed  4> »  ip  given  the 
coefficients  in  the  mapping  function  and  proceeded  to  evaluate  the 
stresses  according  to  equations  (A-4)  and  (A-6). 
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APPENDIX  B 


EVALUATION  OF  STRESS  SYSTEM  1 


Substituting  the  definition  for  4> (p )  in  equations  (9),  inter¬ 
changing  the  order  of  integration,  and  noting  that  e1pZ  = 
e  p<y  [cos  (px)  +  i  sin  (px)]  yields  the  following  (where  z  =  x  +  i  y 
is  a  nondimensional  variable,  not  a  complex  conjugate): 


f°  •  - 

p  JQ(pu)eipz  dp  du 
o 


ra 


Re 


G(u) 


p2J(,(pu)e 


ipz 


dp  du 


r 

p2J0(Pu)e1pZ"  dp  du 

•  0 


(B-l ) 


The  inner  integrals  can  be  evaluated  by  noting  the  following  identities 


r°° 


pJ  (pu)e1pz  dp  =  -  i  5= 
°  dz 


.  d 


r 


J0(pu)e1pz  dp  =  -  i  ~  (uz  -  zz)  2 


p2J0(pu)e"1pZ  dp  =  d 


'o 

2  f°° 


,•  t„ 2  _  —2> 

dz 


(£2 


(B-2) 

J0(pu)e1p7  dp  =  -  (u2  -  z2)  2" 
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where  the  second  integral  has  been  evaluated  using  the  following 
result  from  a  table  of  Laplace  transforms. 


e"st  Jn(at)  dt  =  1 


Substituting  (B-2)  in  (B-l)  yields  the  following: 


gx  +  qv 


ra 


=  Re  i 


dz 


G(u)  du 


gx  _  -  n.  d?_ 


/■a 


—  =yR6^ 


G(u)  du 


—  —  d^ 

Txy =  y Im  ^2 


ra 


G(u)  du 


Equations  (B-4)  can  be  reduced  to  the  well  known  Westergaard 
tion  [27]  for  which  the  stresses  are  given  by: 

ax  =  Re  Z  -  y  Im  Z1 
Oy  =  ReZ  +  yImZt 
Txy  =  -  y  Re  Z' 

by  defining  Z  in  the  following  manner: 

C<T 

.  d  G(u)  du 


(B-3) 


(B-4) 


formul a- 


(B-5) 


(B-6) 
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Substituting  equation  (11)  with  N  =  3  for  G(t)  in  (B-6) . and  carrying 
out  the  necessary  integrations  yields  the  following  expression  for  Z. 

Z  =  -j==5[  +  V'z2'*2  -  a,  -  2Aj  i2  -  f  A5  iA  (B-7) 

where  a-j  =  A-j  +  A3  (2z^  +  a^)/3  +  Ag  (3a^  +  4a2z '  +  8i^)/15 
a2  =  (20A3  z  +  8A5  a2z  +  32A5  z3)/15 

For  the  limiting  case  of  a  through  crack  in  an  infinite  body,  A-j  =  1 , 
A2  =  A3  =  0  and  (B-7)  reduces  to 


which  agrees  with  the  Westergaard  solution  [27]. 


24 


REFERENCES 


[1]  Irwin,  G.  R.,  Discussion  of  the  Paper,  "The  Dynamic  Stress  Distribution 
Surrounding  a  Running  Crack,"  by  A.  A.  Wells  and  D.  Post,  Proceedings, 
Society  for  Experimental  Stress  Analysis,  16,  1,  pp.  69-92  (1958). 

[2]  Fessler,  H.  and  Mansell,  D.  0.,  "Photoelastic  Study  of  Stresses  Near 
Cracks  in  Thick  Plates,"  Journal  of  Mechanical  Engineering  Science, 

4,  3,  pp.  213-225  (1962).  - - 

[3]  Marl  off ,  R.  H.,  Leven,  M.  M.,  Johnson,  R.  L.  and  Ringler,  T.  N.,  "Photo¬ 

elastic  Determination  of  Stress  Intensity  Factors,"  Experimental  Mechanics, 
11,  12,  pp.  529-539  (Dec.  1971).  - 

[4]  Bradley,  W.  B.  and  Kobayashi ,  A.  S.,  "An  Investigation  of  Propagating 
Cracks  by  Dynamic  Photoelasticity,"  J.  of  Experimental  Mechanics,  10, 

3,  pp.  106-113,  (March  1970).  - - 

[5]  Bradley,  W.  B.  and  Kobayashi,  A.  S. ,  "Fracture  Dynamics  -  A  Photoelastic 
Investigation,"  J.  of  Engineering  Fracture  Mechanics.  3.  3.  dd.  317-332 

fr\^  j.  n  r\~n  \  1  '  — — —  **  111  - -  *  r  ~ 


[6]  Kobayashi,  A.  S.  and  Wade,  B.  G. ,  "Crack  Propagation  and  Arrest  in  Im¬ 
pacted  Plates,"  TR-14,  Department  of  Mechanical  Engineering,  University 
of  Washington,  Seattle,  Wash.  (July  1972). 

[7]  Liebowitz,  H. ,  Vanderveldt,  H. ,  and  Sanford,  R.J.,  "Stress  Concentrations 
Due  to  Sharp  Notches,"  Experimental  Mechanics,  7,  pp.  513-517  (1967). 

[8]  Smith  D.  G.  and  Smith,  C.  W. ,  "Photoelastic  Determination  of  Mixed  Mode 

Stress  Intensity  Factors,"  J.  of  Engineering  Fracture  Mechanics,  4,  2, 
pp.  357-366  (June  1972).  - - 

[9]  Smith,  D.  G.  and  Smith,  C.  W. ,  "A  Photoelastic  Evaluation  of  the  Influ¬ 
ence  of  Closure  and  Other  Effects  Upon  the  Local  Bending  Stresses  in 
Cracked  Plates,"  International  Journal  of  Fracture  Mechanics.  6,  3, 

pp.  305-318  (Sept.  1970). 

[10]  Marrs,  G.  R.  and  Smith,  C.  W. ,  "A  Study  of  Local  Stresses  Near  Surface 
Flaws  in  Bending  Fields,"  Stress  Analysis  and  Growth  of  Cracks,  ASTM  STP 
513,  pp.  22-36  (Oct.  1972)~ 

[11]  Smith,  D.  G.  and  Smith,  C.  W. ,  "Influence  of  Precatastrophic  Extension 
and  Other  Effects  on  Local  Stresses  in  Cracked  Plates  Under  Bending 
Fields,"  Experimental  Mechanics.  11,  9,  pp.  394-401  (Sept.  1971). 

[12]  Schroedl ,  M.  A.,  McGowan,  J.  J.  and  Smith,  C.  W. ,  "An  Assessment  of 

Factors  Influencing  Data  Obtained  by  the  Photoelastic  Stress  Freezing 
Technique  for  Stress  Fields  Near  Crack  Tips,"  J.  of  Engineering  Fracture 
Mechanics,  4,  4,  pp.  801-809  (Dec.  1972).  "  ' 


25 


[13] 

[14] 

[1  5] 

[16] 

[17] 

[18] 

[19] 

[20] 


Schroedl,  M.  A.  and  Smith,  C.  W.,  "Local  Stresses  Near  Deep  Surface 
Flaws  Under  Cylindrical  Bending  Fields,"  VPI-E-72-9,  Progress  in  Flaw 
Growth  and  Fracture  Toughness  Testing,  ASTM  SIP  536,  pp.  45-63 
(Oct «  1 973 ) « 

Schroedl,  M.  A.,  McGowan,  J.  J.  and  Smith,  C.  W.,  "Determination  of 
Stress  Intensity  Factors  From  Photoelastic  Data  with  Application  to 
Surface  Raw  Problems,"  VPI-E-73-1,  (Feb.  1  973),  (In  Press)  0.  of 
Experimental  Mechanics.  - 

Harms,  A.  E.  and  Smith,  C.  W.,  "Stress  Intensity  Factors  in  Long  Deep 
Surface  Flaws  in  Plates  Under  Extensional  Fields,"  VPI-E-73-6  (Feb.  1973), 
NCh  Meetin9  °f  Society  for  Engineering  Science,  Raleigh, 

Schroedl,  M.  A., and  Smith,  C.  W.,  "Influence  of  Three  Dimensional  Effects 
Intensity  Factor  for  Compact  Tension  Specimens", 

ASTM ^S TP Pa9es>  (April  1973),  (In  Press)  Fracture 


[21] 


[22] 


Schroedl,  M.  A.,  McGowan,  J.  J.  and  Smith,  C.  W.,  'Use  of  a  Taylor 

upT1?S-n°JIe5f!S!\Mefh0^  In  Photo  elastic  Stress  I  ntens  i  1y  Determinations, " 
VPI -E-73-34  (1  973),  Spring  Meeting  SESA  Detroit,  Mich.  (May  1974). 

Mull  ini  x,  B.  R.  and  Smith,  C.  W.,  "Distribution  of  Local  Stresses 
Across  the  Thickness  of  Cracked  Plates  Under  Bending  Fields,"  (In 
Press)  International  Journal  of  Fracture  Mechanics  (19741. 

McGowan,  J.  J.  and  Smith,  C.  W.,  "Stress  Intensity  Factors  for  Deep 
Cracks  Emanating  from  the  Corner  Formed  by  a  Hole  Intersecting  a 
Plate  Surface, "  VPI -E-74-1  (1  974). 

Smith,  C.  W.,  'Use  of  Three  Dimensional  Photo  elasticity  in  Fracture 
Mechanics,  (Invited  Paper)  ThirVJ  International  Congress  on  Experi¬ 
mental  Mechanics,  Los  Angeles,  Calif.,  May  13-18,  1  973  (In  Press), 

pp^wI^BM^Dec^lDfs)31101  J‘  of  Experimental  Mechanics,  13,  12, 

McGowan,  J.  J.  and  Smith,  C.  W.,  "A  Finite  Deformation  Analysis  of 
VPNE-74  ioel(19747°Und1n9  the  TiP  °f  Crack’Like  Ellipses," 

Savin,  C.  N.,  Stress  Concentrations  Around  Holes.  Pergamon  Press, 
pp.  9-14  ( 1 974 “  - -  8 


[23]  Muskhel ishvili,  M.  I.,  Some  Basic  Problems  of  the  Mathematical  Theory 

of  Elasticity,  Noordhoff,  (1953).  - - - JL 

[24]  Gross,  B.  and  Mendelson,  A.,  "Plane  Elasto-Static  Analysis  of  V- 

Notched  Plates,"  International  Journal  of  Fracture  Mechanics,  8,  3. 
pp.  267-276  (Sept.  1972).  - “ 


26 


[25]  Sneddon,  I.  N.  and  Srivastav,  R.  P.,  "The  Stress  Field  in  the  Vici¬ 
nity  of  a  Griffith  Crack  in  a  Strip  of  Finite  Width,"  Irtt.  J.  Eng. 
Sci.,  9,  pp.  279-488  (1971). 

[26]  Squire,  W..  Integration  for  Engineers  and  ScientlciSts,  American 
Elsevier  Publishing  Company,  Inc.,  New  York,  p.  279  (1970). 

[27]  Westergaard,  H.  M.,  "Bearing  Pressures  on  Cracks,"  Trans.  ASME  (1  939). 

[28]  Sokol nikoff,  I.  S.,  Mathematical  Theory  of  Elasticity,  Second  Edition, 
McGraw-Hill  Book  Co.,  New  York  (1956). 


27 


28 


NORMALIZED  APPARENT  SIF  (KAP/<rmr) 


Figure  3.  Results  of  Elliptical  Notch  Comparison. 
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Figure  4.  Geometries  to  be  Used  in  Notch  Crack  Comparison. 
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NOR 


SQ.  RT  OF  NORM.  DISTANCE  FROM  CRACK  TIP  (r/a) 


Figure  5.  Results  of  Rectangular  Notch  Comparison. 


NORMALIZED  APPARENT  SIF  (KAP/<?1 


SQ.  RT  OF  NORM.  DISTANCE  FROM  CRACK  TIP  (r/a)1 

Figure  6.  Effect  of  Number  of  Terms  on  Rectangular  Notch  Solution 


NORMALIZED  APPARENT  SIF  (KAP/ffj7fa) 


SQ.  RT  OF  NORM.  DISTANCE  FROM  CRACK  TIP  (r/a)'/2 

Figure  7.  Results  of  Rhombus  Comparison. 
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NORMALIZED  APPARENT  SIF  (KAP/?T> 


SQ.  RT  OF  NORM.  DISTANCE  FROM  CRACK  TIP  (r/a)1 


Figure  8.  Results  of  Blunted  Triangular  Notch  Comparison. 


3.0 


v 

Figure  TO.  Solution  for  G(t)  for  Through  Crack  in  Finite  Width  Strip. 
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NORMALIZED  APPARENT  SIF  (KAP/c? 


SQ.  RT  OF  NORM.  DISTANCE  FROM  CRACK  TIP  (r/a)'/2 

Figure  11.  Normalized  Kftp  Curve  for  Through  Crack  in  Finite  Width  Strip 
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